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This paper researches how the systematic errors in phase inversions affect the success rate and 
the number of iterations in optimized quantum random-walk search algorithm. Through geometric 
description of this algorithm, the model of the algorithm with phase errors is established and the 
relationship between the success rate of the algorithm, the database size, the number of iterations 
and the phase error is depicted. For a given sized database, we give both the maximum success 
rate of the algorithm and the required number of iterations when the algorithm is in the presence 
of phase errors. Through analysis and numerical simulations, it shows that optimized quantum 
random-walk search algorithm is more robust than Grover’s algorithm. 
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I. INTRODUCTION 


Quantum computation follows the laws of quantum 
mechanics such as quantum interference, superposition 
and entanglement to achieve a new computing pattern 
of information processing with a powerful parallel com¬ 
puting capability. In 1996, the Grover algorithm^ pro¬ 
posed by Grover is one of the most typical quantum 
algorithms. It only needs 0(y/ 2”) computational com¬ 
plexity to complete a 2 n -sized database search, achiev¬ 
ing quadratic speed-up compared with the classical ex¬ 
haustive algorithm. Grover’s algorithm makes people 
see the great power of quantum algorithms. Since 
then further research based on Grover’s algorithm is 
proposed continuously. I 2 - 3 - 4 ! Grover’s algorithm has also 
been proved to be physically implemented J 5 1 

In practice, errors in the implementation of phase in¬ 
version are inevitable, including systematic errors and 
random errors. The error of each step may be very small, 
but the scale of quantum computation is usually expo¬ 
nential level. So, it may lead to big errors finally and 
affect the algorithm results. The robustness of the algo¬ 
rithm to the error refers to the stability when it is affected 
by error interference. It can be measured by the change 
of success rate before and after the introduction of errors. 
In 2000, Gui Lu Long found systematic errors in phase in¬ 
versions and random errors in Hadmard-Walsh transfor¬ 
mations are the dominant effects in Grover’s algorithmJ 6 1 
In 2003, Shenvi studied the robustness of Grover’s algo¬ 
rithm to a random phase error in the oracle and the error 
tolerability of the algorithm. ! 7 1 It has a significant value in 
studying the robustness of quantum algorithms to errors, 
which can improve the reliability of quantum algorithms 
and benefiting the implementation of algorithms. 


‘Electronic address: 2010thzz@sina.com 


In recent years, quantum random walk has become a 
research hot topic. Quantum random walk is the counter¬ 
part of the classical random walk. It evolves with quan¬ 
tum state and is affected by the quantum interference. 
In 2003, Shenvi proposed a quantum random- walk s earch 
algorithm (SKW algorithm).!®1 It iterates jV2 n+1 times 
and the success rate is 0.5, which is the first algorithm 
using quantum random walk to reach quadratic speed-up 
on an unstructured 2 n -sized database. In 2008, Potocek 
put forward optimized quantum random-walk search al¬ 
gorithm (optimized SKW algorithm) ! 9 1 on the basis of 
the former. It iterates \\[2P times and the success rate 
is 1 approximately, achieving the same calculation re¬ 
sults of Grover’s algorithm. Since then, some quantum 
algorithms based on quantum random walks have been 
proposed. I 10 - 11 - 12 - 13 ! 

In 2006, Li Yun studied the effects of errors in phase 
inversions on SKW algorithm] 14 ! and found that in addi¬ 
tion to reducing the probability of marked state eventu¬ 
ally, the error also reduces the maximum probability of 
the unmarked state, which is different from Grover’s algo¬ 
rithm. While in the relevant experiments (such as NMR), 
as long as the probability of the marked state is much 
higher than other states, the algorithm is still successful. 
One of the main differences between SKW algorithm and 
Grover’s algorithm is the use of the Grover diffusion op¬ 
erator, G = 2 jS’ 0 ') ( S c \ — I. The Hilbert space of SKW 
algorithm is the tensor product H = H c <g> H s , where 
H c is the n dimensional space associated with the quan¬ 
tum coin, and H s is the 2 n dimensional Hilbert space 
representing the vertex. The Grover diffusion operator 

as the coin operator only acts on the n dimensional coin 

n 

space, thus, I*? 0 ) = —k M)- But in Grover’s algo- 
vn d= 1 

rithm, it acts on the 2 n dimensional database, corre¬ 
sponding to the vertex space in SKW algorithm, thus, 
2 n 

\ s °) = £ |a;). The increment of work space will 
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increase the difficulty of quantum gates implementation 
and design costs. Also, it may result in introducing more 
errors, causing more interferences to the algorithm. From 
a practical point of view, the search algorithm based on 
random walk may have more advantages over Grover’s 
algorithm. Therefore, its in-depth study has important 
theoretical significance and practical value. Currently, 
there is no research about the effects of errors on opti¬ 
mized SKW algorithm. 

This paper studies the effects of systematic phase er¬ 
rors produced by the Grover diffusion operator in opti¬ 
mized SKW algorithm. With theoretical analysis and 
numerical simulations methods, we establish the model 
of the algorithm with phase errors and depict the rela¬ 
tionship between the success rate of the algorithm, the 
database size, the phase error and the number of itera¬ 
tions. For a given sized database, the relationship be¬ 
tween the maximum success rate of the algorithm and 
the phase error is given. The relationship between the 
required number of iterations and the phase error is also 
given. A critical value of the database scale is obtained. 
The algorithm will be almost unaffected by the error if 
the database scale is smaller than that. Numerical sim¬ 
ulations show that optimized SKW algorithm is more 
robust than Grover’s algorithm. 

The rest of this paper is organized as follows. In Sec¬ 
tion 2, we introduce optimized SKW algorithm. The ge¬ 
ometric description of optimized SKW algorithm is pro¬ 
posed in Section 3. Section 4 establishes the model of the 
algorithm with phase errors. In Section 5, we give the 
results of numerical simulation and make a robustness 
comparison with Grover’s algorithm. The conclusion is 
summarized in Section 6. 


II. OPTIMIZED SKW ALGORITHM 

Optimized quantum random walk search is based on 
a quantum random walk on a hypercube. The n di¬ 
mensional hypercube is a graph with 2 n vertices, each 
of which can be labelled by an n-bit binary string. Two 
vertices are connected by an edge if their Hamming dis¬ 
tance is 1. The data in the 2 n -sized database is expressed 
in binary, x = (x\, X 2 ■ ■ • x n ), Xi = {0,1}. |x| is its 

Hamming weight and p(x) = \x\ mod 2 is its parity bit. 
Then they can be mapped to a 2 n+1 dimensional space, 
x —> x! = x\\p(x), i.e. adding the parity bit at last. As 
the oracle contains information of the marked state and 
it is one-to-one with x ', we suppose the oracle is still able 
to identify the marked state after mapping. 

When x! corresponds to the vertex of a dimensional 
n + 1 hypercube, the Hamming weights of these vertices 
are even. Assuming that the marked vertex is 10), its 
selection does not affect the final result since the hyper¬ 
cube is symmetrical. This database search problem is 
then equivalent to searching for a single marked vertex 
amongst the even vertices on the hypercube. 

In detail, optimized SKW algorithm can be described 


by the repeated application of a unitary evolution oper¬ 
ator UU' on the Hilbert space H = H c ® H s of n + 1 di¬ 
mensional hypercube. Each state in H can be described 
as \d, x), where d specifies the state of the coin and x 
specifies the position on the hypercube. The operators 
can be written as 


U = S-C 0 , 

(1) 

U' = S-C, 

(2) 

where the shift operator is 


2 n+1 -l n+ 1 

S= ^2 ^ \d,x® e d ) (d,x\. 

x=0 d— 1 

(3) 

It depends on the state of the coin to map a state \d,x) 
onto the state \d,x(Bed), where |e*d) = |0 • • ■ 010 ■ ■ • 0) 
is a null vector except for a single 1 entry at the dth 
component. The coin operator can be written as 

C' = C o ® I + (Cl - C 0 ) ® |0) (0|. 

(4) 


It consists of two parts Co and C\. Co acts on unmarked 
vertices, which is generally chosen as Grover operator. 

i.e. 


C 0 =G = 2\S c )(S c \-I, (5) 


where \S C ) 


71+1 

p= + |d) is the equal superposition 

d—1 


over all n + 1 directions. Ci acts on the marked vertex, 
which is generally chosen as —I. 

Optimized SKW algorithm is implemented as follows: 

(1) Initialize the quantum computer to the equal super¬ 
position over all states, \ipo) = 1+) <8> \S S ). By orthog¬ 
onal projection P e , \ipo) is mapped to the basis vectors 
that Hamming weight of vertices are even to obtain the 


initial state 



(2) Apply the evolution operator UU' = S ■ (Co ® I ) 


S ■ C’ to the initial state 


+o (e) ) with § \/2" times. C' is 


used to flip the phase of the marked state, C' jS"- 7 )® |x) = 
(_!)/(*) | S G ) ® |x), when \x) = |0), f(x) = 1; when 


\x) ± |0>, f(x) = 0. 

(3) Measure the quantum state. The probability of 
marked state) S c ) <g> |0) is 1 — 0(^A_). 

The coin operator C' changes the phase to mark state 
in the algorithm, taking on the function of an oracle. 
It iterates \\[2P times under ideal conditions and the 
success rate is 1 approximately. Each iteration is equiv¬ 
alent to twice walks on the hypercube, containing four 
operations with two coin-flip operations and two trans¬ 
fer operations. C' is only applied once, so each iteration 
only calls oracle once. 

In order to establish and analyze the model of opti¬ 
mized SKW algorithm with phase errors, the geometric 
description of optimized SKW algorithm will be given in 
the following. 
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III. GEOMETRIC DESCRIPTION OF 
OPTIMIZED SKW ALGORITHM 


In order to simplify the problem, let us first show that 
the random walk on the hypercube can be collapsed onto 
a random walk on the line. We use the method presented 
Let Pjj 

Theorem 1 of [8]. 

Theorem 1 U and U' are the unitary operators in 
Section 2 and Pij is the permutation operator. Then 
UU' and Pij satisfy PPUU'Pij = UU'. 


in [8]. Let Pij be the permutation operator defined in 


IJ OO-LIOLJ ij l 3 

Proof : It can be obtained from Theorem 1 in [8] that 
P^jUPij = U and PPU'Pij = U'. We multiply these 
two equations by the left and right respectively and the 
theorem is proved. 

Since the initial state 
Pij with eigenvalue 1 for all i and j. 


*Po 


( e )' 


^ is an eigenvector of 
According to 

Theorem 1, any intermediate state ( UU') 1 ipo' e ^ must 
also be an eigenvector of eigenvalue 1 with respect to 
Pij. Thus (UU'Y preserves the symmetry of ipo^^j 

with respect to bit swaps. Define 2 n basis states, 
\R, 0), \L, 1), |I?, 1) • • • |I?, n — 1}, \L, n), that their spe¬ 
cific forms are the same with equation (10) and (11) in 
[8]. Based on the above basis states, the coin operator 
and shift operator will be rewrote respectively by equa¬ 
tion (13) and (12) in [8]. So, the unitary operators also 
become the form of (14) and (15) in [8]. 

As both U and U 1 are real matrices, UU and UU' are 
also real matrices. Their eigenvalues and eigenvectors are 
complex conjugate pairs. Before analyzing the eigenvalue 
spectrum of UU', let’s analyze the eigenvalue spectrum 
of UU first. 

The eigenvalue spectrum of U is already given in [15]. 
Since UU jug) = U(e ±1Uk jug)) = e ±l2uk juj), we obtain 
the eigenvalue spectrum of UU, 


e ±i2u> k _ cog 2u k ± j gi n 2 oj k 

= 1 + ± Vk(n - k ). (6) 

Theorem 2 U and U' are the unitary operators in 
Section 2. Then there are two eigenvectors of UU' with 
eigenvalues close to 1. 

proof: Supposing the size of the database is 2 n ~ 1 and 
n is a multiple of 4. We construct two approximate eigen¬ 
vectors of UU', ipo^/ and \ipi)> 


2 1 . 
S ( \ 

X=1 


Wi) = l c (Z ■ 

x—0 


n — 1 
\J V 2x 


n — 1 
\ \ 2x + 1 


\R,2x) 


\L, 2x + 2)). 


(7) 


where c = 


\ 


i- 1 

E 

x=0 


n — 1 
x 


. We evaluate 


UU' 




2\/^~T 


obtaining 


tpo 


(e) 


UU' 


^o (e) ) = l- 


( 8 ) 


(9) 


tp is an eigenvector of 


So, when n is large enough, 

UU' with eigenvalue close to 1. Then we evaluate 

= (|l?,f-2) + |L,f +2» 


UU'\ipi) 

= i^)-|(i -l) 1 


12 
c n 


\ 


n — 1 

N V f -1 


n— 1 

n _ 2 

(|7f> + |7f» 


in the same way to obtain 

{ipi\UU'\'ip 1 ) = 1- 


1 


n — 1 
72 


( 10 ) 

( 11 ) 


Since 1 < c 2 < 1 H—, |^i) is also an eigenvector of UU' 
with eigenvalue close to 1. The theorem is proved. 

Then according to Theorem 2,3 in [8] and the eigen¬ 
value spectrum of UU, it can be obtained that there are 
exactly two eigenvalues of UU' with their 

real parts bigger than 


L = 1 - 


8 n 


3(n + l) 


2 ' 


( 12 ) 


Set the corresponding eigenvectors |w'o) and |— ui'o), sat¬ 
isfying |u/o) = |— oj'o)* . According to the Theorem 4 in 
[8], |w'o) an d |— uj'o) can be also proved to be approx¬ 


imately composed of the linear combination of 


ipo 


(e) 


f n- 1 \ 



X '' \R,2x) 




ipo 


(e) 


* l^l)), 



( n-1 \ 


+\ 

- -\L, 2x)), 



(13) 


( 14 ) 
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Finally, we evaluate the range of u/ 0 . First we evaluate 
to obtain 


(t/>i|C/[/'|^o (e) ) = (Vi 


i>o 


(e) 


+ 


1 2 


V’o 


(e) 


uu' I^i) = (v>o (e) | ih)~ 


c x/2"’ 
1 2 


a/2"’ 


(15) 


(16) 


Supposing |a) = ^=( + i|i/d)), according to the 

Theorem 5 in [8], 

\e iuj '° - (a| U' |a) | 


< 2(^1 
— v n* 


3(n+l) 2 


2"+ 4 


3(n+l) 2 


)=0(1L^) (17) 


16c 2 -n- 


2 + 1 
2 ' 


is obtained. The imaginary part is 
Im(a|t/» = i((^o (e) UU' \ip\) — (V'i| UU' V>o (e) )) 


2 

c V2" ' 


On the other hand, because of 

| sinwQ — Im(a| U' |a) | = |Im(e*“ 0 — (a| U' |a))| 


(18) 


3 /2 

,n . 


<\e i “°-(a\U'\a)\<0(—) (19) 


SO, 


, • / 2 , n^2 

sinwj H-— < 0(——). 

0 c- \/2" “ 2" ; 


( 20 ) 


When Wq tends to 0, there will be sinwg = u>' 0 + O(lu' 0 ). 
So 


2 „.n 3/ 2 , 2 

-7= - Of——) <w <-= 

c . V 2 n ’ ~ 0 - c • y/2" 


n 3 /2 

+ 0(^r). (21) 


At last this algorithm can be geometrically described 
through the above conclusions. Starting with the initial 


state 


</>o (e) ), it successively iterates t times, 


{UU') 1 |^o (e) ) = ^(e ia/ot |u/ 0 > + e~ iu ' ot |-w'o» 


= COSWgt 


V>o (e) ) - sinwgf l^i). (22) 


Each iteration is equivalent to a rotation of angle |ojq | 

( e ) \ _ 


on the plane composed of the initial state V'o 
7 f(l w, o) + |-w , o)) and |Vh) = ^(|w'o) - |-w' 0 )) (shown 
in Fig. 1). Because of (ip= 0, the required 

iteration times is 2 \Z'o\ = f V^- After obtaining \ipi), 
the marked state | R, 0) can be measured with probabil¬ 
ity p = |(i?,0 | = 2^. So the success rate of the 

algorithm is 1 — 0(+j-). Both the number of iterations 
and the success rate meet the conclusions in [9]. The 
geometric description of optimized SKW algorithm has 
been presented. 



FIG. 1: Each iteration of the algorithm UU' can be viewed as 
a rotation on the plane composed of the initial state |-0c/ e ^ 
and the final state \ipi). The rotation angle is |wo| 


IV. MODEL OF THE ALGORITHM WITH 
PHASE ERRORS 

A. Build the model 

The model of the algorithm with phase errors can be 
built after the geometric description of the original al¬ 
gorithm. Let the database size be 2 n ~ 1 . As the Grover 
operator imperfections in phase inversions of the algo¬ 
rithm are systematic, the explicit form of the inaccurate 
operator can be written as 

Co = G={l-e ie )\S c )(S c \-I, (23) 

where 6 = n + S and S is a constant error. When 5 = 0, 
it recovers the original operator. Then the coin operator 
with errors can be written as 

C' = Co ® I + (Ci - Co) ® |0) (0|. (24) 

The evolution operator becomes 

U = S-C 0 , (25) 


U' = S-C'. (26) 

Each iteration of the algorithm becomes UU'. 

We analyze the eigenvalue spectrum of UU first. The 
eigenvalue spectrum of U is already given by [14] with 


e iUk = l(l-e ie )-(l-e i0 )-T^ 

l n 2 


'a 


■ +4e ie , 
(27) 
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k 7 0, k When k = 0, e lu, ° = —e' e = cos 5 + 7 sin (5; 
k = n, e laJn = e“. The introduction of errors lead 
to an <5 angle deviation between e luJ ° and e*“° = 1 
on the polar coordinate plane (see fig. 4.1 in [17]). 
The eigenvalue spectrum of UU is e 2luJk , when k = 0, 
e 2 “° = e 2 ' 6 = cos 26 + i sin 28. So, there is an 2<5 angle 
deviation between e 2lw ° and e 2lU10 on the polar coordi¬ 
nate plane. Eigenvalues are no longer complex conjugate 
pairs. 

For the eigenvalue of UU', it only needs to analyze 
e l “ ,Q , e IU ' 1 I which originally correspond to the eigenvec¬ 
tors |w' 0 > and |— ui’o) , because they determine the initial 
state and the final state of the algorithm. These two 
eigenvalues are conjugate pairs originally. But the intro¬ 
duction of errors lead to them being no longer conjugate. 
When 8 increases, the difference of these two eigenvalues 
gradually increased. e lU1 ' 0 and e lu ' 1 will gradually ap¬ 
proach 1 and e l2e respectively according to [14]. If <5 is 
very small, they are still in the vicinity of 1 . 

According to Theorem 1 and [14], the random walk 
which introduces phase errors can also be collapsed onto 
a random walk on the line. To redefine the operator, 
where 

U' = U-(l-e ie )\L,l)(R,0\, (28) 

the specific form of U has been redefined in [14]. 

Expanding eigenvalues of UU to the second order with 
respect to 8 and using the Theorem 2 and Theorem 3 in 
[ 8 ], it is can be proved that there are exactly two eigen¬ 
values of UU' with their real parts L bigger than 


1 - 



4)_ (2- i): 


3(1 +Jr 



(29) 

When (5 = 0 and n is replaced by n + 1, it is equal to 
equation (13). Then we evaluate 


so 


<V> 1 |t/i/'|V’i) = i- 


(e ie - l)e 


2 c 2 


71—1 

72 


= 1 - 


(e iS + l)e i( 


2 c 2 


71—1 

72 


(33) 

From the above equations (32) and (34), when n is 


big enough, 


is an eigenvector of UU' with eigen¬ 
value close to e 2lS and |i/q) is an eigenvector of |i/ , i) with 
eigenvalue close to 1. Because e 2lS = cos2(5+isin2(5, its 
real part can be expanded by Taylor expansion cos 2<5 = 
1 — 2(5 2 . Based on the analysis above, when the error is 
small, the algorithm can still be viewed as rotation on 
the two-dimensional plane composed of (7) and ( 8 ). 


B. analysis of the algorithm 


With the error increasing gradually, uj'q j and w'i 

are no longer conjugate. They can’t be expressed as 
the linear combination of (14) and (15). So the algo- 

if >(77 to \ipi) completely, and 


rithm can’t rotate from 
the probability to obtain |'i/q) becomes small finally. 


Expanding 
spectively, we obtain 


tpo ^ ) and with 


<*/ o ) and 


u>' i) re- 


'0o (e) ^ = a 0 u/o^+ai W, i) + \A a o ~~ a i l7o), (34) 


l^i) = b 0 


w' 0 


w, i) + V 1 “ 6 o~ 6 il7i), (35) 


where I 70 ) and I 71 ) are the normalized vectors orthog- 


oj'o) and 


onal to 


numbers, decided by S and n. Evaluating 


u>' 1 ). ao,ai,bo,bi are the complex 


UU' 


(e) \ _ e 2i0 


777 =e 
+ 


*> w >-:*&( 


1-e™ (l-e™-n 


V2 n ~ 1v n 


170) 


(l-e ig )-x/n^T 


+7)), 


so 


V’O 


(e) 


UU' 


*,<■>)= 

__ 2i5 _ (e^ + l)e^ 

— e On j 


(30) 

(31) 


In the same way 


1 - 5 ^ ~ <5 2 (2 - Jr) 

= Re((Vo 


» 


= cos <- 


UU' 


W e) )) 

7o (e) ' 2 


= Oq cos wo + a 2 cos u>\ + ]+ cose 
~ i#o,i 

< Oq + a\ + L( 1 — Oq — al) , 


we obtain 




4>o 


(e) 


(36) 


UU' \ipi) 

= hM-^(i-*)- 


1 / 1—e 


20 (l—e id ) 2 \ 


n — 1 

N V f - 2 


a 2 + a 2 >l--^ T [-l I + (2-| r )(5 2 ]. (37) 

(|7f-2) + |L,f+2)) 1 -L 2 2 

In the same way, evaluating Re (ipi\ UU' \ipi), we obtain 


c ' 4 


2 n 


n — 1 

NV1-1 


(|7f) + |7f)), 


,2 , 7 2 . n 1 / 4 5(52 (38) 


( 32 ) 


b 2 +bl> l--7,( 

1 - 1j 


4 c 2 


71—1 

72 
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Now the algorithm in the presence of phase errors can 
be geometrically described through the above conclu¬ 


sions. Applying the operator on the initial state 
successively with t times, we get 


V’o 


(e) 


( UU'Y 


V’O 


(e) 


= a 0 e 

_ aobie ZUJ Qt — aiboe luJ 1<: 

aobi —aibo 

CLQa\(e XU} o t _ e lCt ’ 
aobi— aibo 


w'o) + aie laj,lt w\) + yjl - al~ a\ \^ 0 ) 


^o (e) )-^l7o)) 

(IV’i) - y/ei 111)) + V^o iTo) > 

_ _ (39) 

w here = (C/C/')* Ito>, a / 1 “ a o “ a i = \/^o and 
\/l — — = VeT- So, the amplitude of \ipi) is 

tn = |S (40) 

a 0 Oi — ai6o 


decided by eigenvalues e lw 0 , e luJ 1 , iteration times t and 
coefficients ao,ai,&Oj&i- 

When 5=0, there are e lu '° = e lw °, e lw>1 = e _J “°, 
a 0 =ai=4|, b 0 = - bi = A= and e 0 ,£i ->-0. So w = 


sin wq t, it means that 


can completely rotate to 
|^i) finally. The probability to obtain \ip\) is nearly 1. 
When 5 ^ 0, £o,£i are no longer close to 0 and there 
is another state other than \ipi) at last. So the proba¬ 
bility to obtain \ipi) will decrease. On the other hand, 


(e) 


(Vi V’o (e) ) = 0, ipo 

the plane composed of 


on 


needs to rotate angle of f 
^o (e) ) and \ipi) in theory. It is 


equivalent to making uj'ot — u'\t 


= 7r at last. When S 


increases, the difference between eigenvalues ui'q—uj'i will 
increase. So, the number of iterations t to reach maxi¬ 
mum w will decrease. There is e lt in the equation (41), so 
it is a periodic function. Therefore when the algorithm is 
in the presence of phase errors, its success rate will still 
periodically change with the number of iterations. 

The above analysis shows that, after the introduction 
of errors, the success rate of the algorithm decreases and 
the required number of iterations decreases too. The suc¬ 
cess rate still periodically changes with the number of 
iterations. 


V. NUMERICAL SIMULATIONS AND 
ROBUSTNESS COMPARISON 

Through the above analysis and proof, the relationship 
between the size of database, the phase error, the num¬ 
ber of iterations and the amplitude of the marked state is 
obtained. We also simply analyze the effects of the phase 
error on the algorithm. In this section, we will respec¬ 
tively analyze the impact of phase errors on the success 
rate and the number of iterations through numerical sim¬ 
ulations. Finally, the quantitative relationship between 


the size of database, the phase error, the number of iter¬ 
ations and the success rate is given. 

Firstly, we simulate the relation between the number of 
iterations and the success rate when n = 8. We compare 
the two cases, where 5 = 0 and 5 = 0.2, as shown in Fig. 
2 . 



FIG. 2: The relation between the success rate of the algorithm 
p and the number of iterations t when n = 8. The solid line 
represents 5 = 0.2 and the dotted line represents 5 = 0 


The solid line indicates that the algorithm without er¬ 
rors iterates 13 times to reach the maximum probability. 
The dotted line represents that it needs 6 iteration times 
to reach the maximum probability when the error is 0.2. 
But this time the maximum probability decreases com¬ 
pared to the former. Both the success rates periodically 
change with the number of iterations. Simulation results 
are consistent with the theoretical analysis in Section 4.2. 

Then we study the relationship between the maximum 
success rate of the algorithm, the size of database and 
the phase error. Each curve represents different phase 
errors, with 5 = 0,5 = 0.01,5 = 0.001 and 5 = 0.0001, as 
shown in Fig. 3. 

It is shown in Fig. 3 that the maximum success rate 
increases with the increment of the database size first. 
When it reaches a threshold, the maximum success rate 
decreases exponentially. And the larger the error is, the 
smaller the database scale will be when the maximum 
success rate begins to decrease. With curve fitting, we 
obtain the expression about the maximum success rate 
Pmax, the size of database 2 n and the phase error 5, 


3.8p 0 

Pmax ~ 3.8 + 5 2 2«- 


(41) 


Po is the success rate of the algorithm without errors. It 
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FIG. 3: The maximum success rate of the algorithm p max vs 
the scale of the database n with different phase errors. The 
asterisks line for 5 = 0, the + line for 8 = 0.0001 , the circles 
line for 8 = 0.001 , the x line for 8 = 0.01 


FIG. 4: The scale of the required iteration times log 2 t vs 
the scale of the database n with different phase errors. The 
asterisks line for 5 = 0, the + line for 5 = 0.0001 , the circles 
line for 5 = 0.001 , the x line for 5 = 0.01 


can be obtained from Section 4 


Po 


1 

n /2-i 

y - 



(42) 


Combining the above two equations (42) and (43), we 
obtain 


Pmax — 


19 

c 2 (19 + 55 2 • 2”)' 


(43) 


It can be known from equation (44) that before the 
search in order to ensure the success rate of the algo¬ 
rithm, we need to predict the size of the phase error to 
limit the size of the database. For example, if the re¬ 
quirements of the success rate is not less than 0.5, the 
database size can’t be more than 

Then we study the relationship between the iteration 
times required to reach the maximum success rate, the 
size of the database and the phase error. The correspond¬ 
ing simulation results are shown in Fig. 4. 

It is shown in Fig. 4 that the required iteration times t 0 
increases with the increment of the database first. Then 
it remains the same. When the required iteration times 
begin to remain the same, the larger the error is, the 
smaller the database scale will be. With curve fitting, we 
obtain the expression about iteration times the size 
of the database 2 n and the phase error 5, 


exceed a critical value, the maximum success rate of the 
algorithm and the required number of iterations are the 
same with the original algorithm respectively. It suggests 
that this algorithm has certain tolerance to this kind of 
errors. Comparing Fig. 3 with Fig. 4 for the same error, 
we find that the corresponding database scale when p max 
begins to decrease is smaller than that when t a begins 
to remain the same. It is shown p max is more sensitive 
to the influence of the error. So, we take the value of n 
when p max begins to decrease as the critical value. p max 
is already obtained. When p max begins to decrease, we 
obtain the corresponding (5, n) by numerical simulations. 
The expression about the database scale n and the phase 
error 8 is 

0.4642 

n = 1.806 • log 2 6 . (45) 

So, when the algorithm is in the presence of errors, it is 
almost unaffected as long as the database scale satisfies 

0.4642 

n < 1.806 • log 2 5 . 

Now we respectively obtain the maximum success rate 
Pmax and the required iteration times t 0 expressions 
about the phase error 8 and the database size 2". From 
the analysis in Section 4.2, the probability to obtain 
marked state p is a periodic function, which can also be 
confirmed from Fig. 2. So, making use of equation (44) 
and (45), we obtain an expression of the database size 
2", the phase error 5, the number of iterations t and the 
success rate p, 


t 


O 



(44) 


19 


P = 


c 2 (19 + 55 2 • 2 n ) 


sin (\/™ + 62 •*)• ( 46 ) 


From the above analysis, when it is in the presence 
of phase errors, as long as the database scale does not 


According to equation (47), before the search it can 
depend on the requirement of the success rate to weigh 


























the number of iterations according to the given sized 
database and the predicted error. It will improve the 
algorithm efficiency. It is shown in Fig. 5 when n = 8. 
When <5 = 0 or <5 = 0.2, it is observed exactly the same 
as that in Fig. 2. 

n=8 




FIG. 6: the probability gap of optimized SKW algorithm Api 
vs the scale of the database n with different phase errors. The 
asterisks line for 5 = 0, the + line for 5 = 0.0001 , the circles 
line for 5 = 0.001 , the x line for S = 0.01 


FIG. 5: The three-dimensional graph of equation (47) when 
n = 8. The x-axis is the number of iterations t . The y-axis 
is the phase error 5 and the z-axis is the success rate of the 
algorithm p 


Finally, we make a simple comparison between opti¬ 
mized SKW algorithm and Grover’s algorithm of robust¬ 
ness to the phase error. The definition of probability gap 
is the probability difference between the marked state 
and the maximum rest state! 16 ! it is defined as 

Ap = p out - max(p 2 ,P3 • • -PN-i)- (47) 

In practice, as long as there is sufficient probability gap, 
the probability of finding the marked state is much higher 
than other states at the end of the algorithm. It is still 
successful. If the algorithmic probability gap varies little 
before and after the presence of errors, its robustness is 
strong. So, we can measure the difference between their 
probability gaps before and after the presence of errors 
to compare robustness. 

According to [17], the expression of the database size, 
the phase error and the maximum success rate is given 
in Grover’s algorithm, 


4 . . 

Pmax “ 4 + §22n ■ ( 48 ) 

So, the probability gap of Grover’s algorithm is 

A B =p„«-T^p. (49) 

Firstly, we simulate the probability gap of optimized 

SKW algorithm Api with different phase errors. It is 
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FIG. 7: the difference between the probability gaps Api — Ap 2 
vs the scale of the database n with different phase errors. The 
asterisks line for <5 = 0, the + line for S = 0.0001 , the circles 
line for S = 0.001 , the x line for 6 = 0.01 


shown in Fig. 6. Then we simulate the difference between 
the probability gaps Api — Ap 2 5 shown in Fig. 7. 

There is a tendency that the probability gap of opti¬ 
mized SKW algorithm decreases exponentially with the 
increment of errors, shown in Fig. 6. But for the same 
size of databases the difference between the probability 
gaps increases gradually, as shown in Fig. 7. This means 
that the probability gap of Grover’s algorithm decreases 
faster, and optimized SKW algorithm is more robust. 
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VI. CONCLUSION 

This paper studies the effects of systematic errors in 
phase inversions on optimized SKW algorithm in detail. 
The geometric description of the algorithm is given out 
first. Then the model of this algorithm with phase er¬ 
rors can be built. It is found that phase errors will 
reduce both the number of iterations and the success 
rate of the algorithm through theoretical analysis and 
numerical simulations. The larger error will lead to the 
smaller maximum success rate of the algorithm and the 
fewer required iteration times. We obtain two expres¬ 
sions about the maximum success rate of the algorithm 
and the required iteration times, which are both related 
to the phase error and the database size. Then we ob¬ 
tain a critical value of the database scale n. When 


0.4642 

n < 1.806 • logj^ , the algorithm is almost unaffected 
by the error. Based on the above relationships, we de¬ 
pict the relationship between the success rate of the al¬ 
gorithm, the phase error, the number of iterations and 
the database size. Finally, Analysis and numerical sim¬ 
ulations show that optimized SKW algorithm is more 
robust than Grover’s algorithm. 
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